Chaotic imaging in frequency downconversion 

Emiliano Puddu, Alessia Allevi and Alessandra Andreoni 
Dipartimento di Fisica e Matematica, 
Universitd degli Studi delVInsubria, via Valleggio, 
11, 22100 Como, Italy and I.N.F.M., 
Unita di Como, via Valleggio, 11, 22100 Como, Italy 

Maria Bondani 

I.N.F.M., Unita di Como, via Valleggio, 11, 22100 Como, Italy 

Abstract 

We analyze and realize the recovery, by means of spatial intensity correlations, of the image 
obtained by a seeded frequency downconversion process in which the seed field is chaotic and an 
intensity modulation is encoded on the pump field. Although the generated field is as chaotic as 
the seed field and does not carry any information about the modulation of the pump, an image of 
the pump can be extracted by measuring the spatial intensity correlations between the generated 
field and one Fourier component of the seed. 
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During the past twenty years, much theoretical and experimental work has been devoted 
to the study of image transfer and non-local imaging. The spatial properties of nonlinear 
interactions have been largely used for the realization of such experiments From 
a quantum point of view, the generation of couples of entangled photons by spontaneous 
downconversion has allowed, through coincidence techniques I2J], to transfer spatial informa- 
tion from one of the twin photons to the other SMB. The same techniques yielded 
similar results with a classical source of correlated single-photon pairs j?l Isj]. It has been 
theoretically shown that, also in the many photon regime, image transfer can be imple- 
mented with both quantum entangled B^™^, mUm. One 
may alternatively place the object to be imaged on the beam pumping the spontaneous 
parametric downconversion process [l^. Images have been actually detected by mapping 
suitable photon coincidences between the generated beams 14, 3|- The corresponding 
imaging configuration in the many photon regime has not been studied yet. In this Letter, 
we propose a classical experiment whose aim is to mimic the many-photon quantum exper- 
iment. We realize a scheme of seeded parametric generation in which an image is encoded 
on the pump field, E 3 (object field), and the seed field, Ei, is chaotic, being obtained from 
a plane wave randomized by a moving diffusing plate. We show that the generated field, 
E2, is also chaotic and does not carry information about the image encoded on the pump. 
This interaction somehow simulates the corresponding spontaneous process, since Ei has 
the same chaotic statistics of the states produced by spontaneous parametric generation. 
Actually, the chaotic field E2 reconstructs an incoherent superposition of images of the ob- 
ject generated by the plane- wave components of Ei 12, Here we show that it is 
possible to recover a single image by measuring a suitable spatial intensity correlation func- 
tion between the seed/reference field Ei and the generated field E2, thus suggesting a way 
to extend quantum imaging protocols to the continuous- variable regime, with the definite 
advantage of measuring intense signals instead of single photons. 

As shown in FigUJa), the amplitude modulation Uo (%o,yo) of field E 3 is obtained by 
placing an object-mask, O, on the beam path. The lens that is located on plane (xl,Vl), 
at distance dp before the nonlinear crystal, images O into O' on plane (xi,yi) at distance 
d = 2f — d F beyond the crystal (2/ — 2f system). A plane- wave seed E 1; slightly non- 
collinear to the pump (not shown in the figure), would generate an E 2 field reconstructing 
a real image of O' on plane (£2,2/2) at a distance S2 = (^2/^3)^ along the E 2 propagation 
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direction, where kj are the wave- vector magnitudes [18[- The amplitude Uf on the plane of 



the crystal entrance face (xf,Vf) is 
U F (xf,Vf) 



2irid 

x / dx dy U (x ,yo) e" 2d > 



For a non depleted pump, the field complex amplitudes 01,2,3 (L) at the crystal exit face are: 



ai(L) = ai(0)cosh(/(7?,^)|a 3 (0)|L)~ai(0) 



a 2 (L) = ig eff Lal(0) 



a 3 (0) 



o 3 (0) I 

ig eff Lm/3)al{0)a 3 {0) 



smh(f(0,P)\a 3 (O)\L) 



(2) 



a 3 (L) = o 3 (0) = C/p (x F ,y F ) , 

where L is the crystal thickness, /($, /3) is a function of the propagation angles (see Fig. [T] 
b)), (7 e // is the coupling constant of the interaction and the approximations hold in the low 
gain regime [19j. 

Now we consider the interaction that occurs with a chaotic Ei (r) cx 
J2n=i a i," ex P ( — * ^i," ' r ) having random complex amplitudes, ai <n , and wave vec- 
tors, ki, n , with random directions but equal amplitudes, k\ )Xl = 2tc/\i, for an ordinarily 
polarized seed wave. Inside the nonlinear medium, each of the spatial Fourier components 
of the seed field that is phase matched with the pump generates an a 2 (L) contribution 
according to Eq. (J2J). Thus N images are simultaneously generated and the overall field is 



N 



E 2 (iw) oc ^^ e// L/(^,/3)at in (0)a 3 (0)e- 
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in which r out = (x out , y ou t, L), see Fig. HJa). Since the wave vectors k 2j „ are linked to ki )n and 
k 3 by the phase- matching condition (k 3 = ki jn + k 2in ), they have random directions, thus 
impairing the reconstructed image hologram visibility. If we let field E 2 propagate freely to 
the plane (x 2 ,y 2 , {k 2 /k 3 )d) where all N images do form for a type-I interaction of paraxial 
beams, we can calculate the intensity on that plane, which turns out to be 18, 3] 

J 2 (r 2 ) cx |E 2 (x 2 ,y 2 , {k 2 /h)d)\ 2 



A' 



X 2 ,n X 2 , y 2) n 
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where the coefficients \c n \ summarize all constant factors and the transverse translations, 
x 2>n and ?/2,n (if refraction at the exit face of the crystal is neglected: x 2 ^ n = (k 2 /k 3 )d sin /3 2 ,n 
and 7/2,n = (^2/^3)^ cos sin 7?2,n) see Fig- d b)) are due to the different directions of the 
k2,n wave vectors. In the last line of Eq. (J3J) the intensity I 2 {x-i) has been written as the sum 
of N terms because, owing to the incoherence of the N components of E 2 , all the interference 
terms vanish. The observation that each term of the sum is proportional to the intensity 
of the n-th reference field component, \ai jn \ , provides a means to recover the image. In 
fact, by evaluating the spatial correlation function of the intensity I\j = |oij| of a single 
component of the seed with the intensity map of the generated field, I 2 (x 2 , y 2 ,z 2 ), we get: 

G (hj, h) = (hjh) - (h^ih) 

N 

2 2 1 2 1 2 

= (Kjl 2^l c «l l a i,nl \Uo(x 2tn -x 2 ,y 2tn -y 2 )\ ) 

n=l 
N 

- (\ a i,j\ 2 )C^2 i c ™i 2 i ai ' n ' 2 \ u ° ( X2 ' n ~ x2 ^ 2 < n ~ ^)i 2 ) 

71=1 

N 

= ^ \ c n\ 2 Wo (>2,n - X 2 ,y 2<n - 7/ 2 )| 2 

71 = 1 

x ((Kjl 2 |ai in | 2 ) - (Kj| 2 )(K,n| 2 )) 

N 

= ^2 l Crl ' 2 \ U ° ( X2 > n ~ X2 '^ 2 ^ _ 2/2) | 2 0" 2 (kl,n| 2 ) $j,n 
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= \cj\ 2 cr 2 (h,j)\U (x 2J - x 2 ,y 2J - y 2 )\ 2 , (5) 
where (• • • ) is the ensemble average. In deriving Eq. (jSJ) we neglected the pump fluctuations 

2 2 i2\/i 2 2 

and used the property of chaotic light: (|ai (J -| |oi, n | ) — (l a i,jl )(l a i,n| ) — ° (l a i,n| )$j,n, 
in which cr 2 (|aij| 2 ) = a 2 (hj) is the variance. This result shows that G I 2 (x 2 ,y 2 )) is 
proportional to the intensity map of the difference-frequency generated image of the object 
\Uo\ "reconstructed" by the j-th component of the seed field Ei. For an experimental 
proof, we have to map I 2 (r 2 ) and to measure I lt j, which is easily done by mapping Ji (ri) 
in the focal plane of a lens. 

The experimental setup is sketched in Fig. El The wavelengths of the interacting fields 
are Ai = A2 = 1064 nm, A3 = 532 nm. Pump and seed fields are obtained from a Nd:YAG 
laser (10 Hz repetition rate, 7-ns pulse duration, Spectra-Physics). The nonlinear crystal 
is a type I /5-BaB 2 04 crystal (cut angle 32°, 10 mm x 10 mm x 4 mm, Fujian Castech 
Crystals). The detection planes of Iij and I 2 (r 2 ) are made to coincide on the sensor of 
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the same CCD camera (Dalsa CA-D1-256T, 16 /im x 16 /zm pixel area, 12 bits resolution, 
operated in progressive scan mode), so that each signal occupies half sensor. The chaotic 
field Ei is generated by passing the beam at Ai through a ground-glass diffusing plate, which 
is moved shot by shot, by selecting a portion of diffused light with an iris, PH, of ~ 8 mm 
diameter and finally by filtering the ordinary polarization with a polarizing beam splitter, 
PBS, and a half- wave plate. Lens L2 (/ = 15 cm) provides the Fourier transform of Ei. To 
check the chaotic nature of this seed field we measured the probability distribution, P r (ii), 
of the intensity recorded by the different CCD pixels for a single shot, and the probability 
distribution, Pt(h), of the intensity recorded by a single pixel for many successive laser 
shots (see Fig. |2l inset). The good agreement with thermal distributions shows that Ei is 
actually randomized in space at each shot and that any Iij is random from shot to shot. 
The desired imaging configuration was realized by using a copper sheet with three holes 
(~ 256 /im diameter) as the mask producing object O, and by recording I 2 (r 2 ) at the 
laser repetition rate. The intensity correlation function in Eq. (JHJ) was evaluated over 1000 
shots by taking the whole map I2 (22,2/2) and by selecting the value of a single pixel in 
the intensity map of I\. In Fig. El a), we show the resulting reconstructed image (map of 
G (h,j, h{x2, 2/2))), to be compared with the plane-wave image obtained in single shot upon 
removing the light diffusing plate, Fig. 01 b). The similarity in the quality of the two images 
is really impressive, in particular if the reconstructed image is compared with a single-shot 
intensity map J 2 (22,2/2), Fig. 01c), and with the average intensity map (I2 (22,2/2)) of the 
1000 repetitions, Fig. Old). 

In conclusion, we have demonstrated that the spatial intensity correlation properties of 
the downconversion process can be used to recover a selected image from a chaotic ensemble. 
The image recovered by G 2/2)) fulfils the properties of the difference-frequency 

generated image that would be obtained by using the single plane- wave Eij as the seed field. 
We expect that the method also works in the case of an unseeded process in the continuous- 
variable regime, in which the selection of a single spatial and temporal frequency in the 
parametric fluorescence cone should determine the position of the reconstructed image. 

The authors thanks A. Gatti (I.N.F.M., Como) for stimulating discussions, 
I.N.F.M. (PRA CLON) and the Italian Ministry for University Research (FIRB n. 
RBAU014CLC_002) for financial support. 
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List of Figure Captions 

Fig. a) Propagation scheme; NLC, nonlinear crystal, b) Interaction inside the crystal, 
optical axis on the shaded plane. 

Fig. |21 Experimental setup: HS, harmonic separator; D, diffusing plate; Mi_5, mirrors. Lens 
Li images O into O' through a 2f — 2f system. Lens Li_2, lenses. Distances are in cm. 
Inset: logarithmic plot of P r;t (Ii). 

Fig. EJ a) Map of G (Ix,j, h{x2, V2)) evaluated on 1000 shots, b) Plane wave image. Chaotic 
images: c) single-shot and d) average over 1000 shots. 
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FIG. 1: a) Propagation scheme; NLC, nonlinear crystal, b) Interaction inside the crystal, optical 
axis on the shaded plane. 
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FIG. 2: Experimental setup: HS, harmonic separator; D, diffusing plate; Mi_5, mirrors. Lens 
Li images O into O' through a 2/ — 2/ system. Lens Li_2, lenses. Distances are in cm. Inset: 
logarithmic plot of P r; t(Ii)- 
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FIG. 3: a) Map of G (Iij, li^2i 2/2)) evaluated on 1000 shots, b) Plane wave image. Chaotic 
images: c) single-shot and d) average over 1000 shots. 
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